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ABSTRACT 



Using wide-field BVRcI imaging for a sample of 16 intermediate redshift 
(0.17 < z < 0.55) galaxy clusters from the Canadian Network for Observa- 
tional Cosmology (CNOCl) Survey, we investigate the dependence of cluster 
galaxy populations and their evolution on environment. Galaxy photometric 
redshifts are estimated using an empirical photometric redshift technique and 
galaxy groups are identified using a modified friends-of-friends algorithm in pho- 
tometric redshift space. We utilize the red galaxy fraction {fred) to infer the 
evolutionary status of galaxies in clusters, using both individual galaxies and 
galaxies in groups. We apply the local galaxy density, S5, derived using the fifth 
nearest-neighbor distance, as a measure of local environment, and the cluster- 
centric radius, tcl-, as a proxy for global cluster environment. Our cluster sample 
exhibits a Butcher-Oemler effect in both luminosity-selected and stellar-mass- 
selected samples. We find that fred depends strongly on S5 and tcl, and the 
Butcher-Oemler effect is observed in all S5 and tcl bins. However, when the 
cluster galaxies are separated into tcl bins, or into group and non-group sub- 
samples, the dependence on local galaxy density becomes much weaker. This 
suggests that the properties of the dark matter halo in which the galaxy resides 
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have a dominant effect on its galaxy population and evolutionary history. We 
find that our data are consistent with the scenario that cluster galaxies situated 
in successively richer groups (i.e., more massive dark matter halos) reach a high 
fred value at earlier redshifts. Associated with this, we observe a clear signature 
of 'pre-processing', in which cluster galaxies belonging to moderately massive 
infalling galaxy groups show a much stronger evolution in f^ed than those classi- 
fied as non-group galaxies, especially at the outskirts of the cluster. This result 
suggests that galaxies in groups infalling into clusters are significant contributors 
to the Butcher-Oemler effect. 

Subject headings: galaxies:clusters:groups:general:evolution 



Introduction 



With increasingly detailed studies of galaxy properties such as stellar population and 
morphology, our knowledge of galaxies in the Universe has greatly improved. The large- 
scale structure in the Universe shows that some galaxies are located in low-density filaments 
while others populate high-density clusters. This complicates our understanding of galaxy 
evolution over the Hubble time, because galaxy properties depend not only on time but also 
on environment. Therefore, in order to understand galaxy evolution in the Universe, galaxy 
properties require viewing simultaneously and systematically over time and in different en- 
vironments. 



The morphology- density relation flDresslerlll980l ) is well- accepted evidence for the envi- 
ronmental dependence of galaxy properties. Early-type galaxies favor dense environments, 
while late-type galaxies are common in less dense environments. Following this pioneering 
work, a considerable number of studies on this topic have been carried out or e xtended (e.g.. 



Dommguez et al.ll2002l : iGoto et al.ll2003l : iHogg et al.ll200a iTanaka et al.ll2005l : iNuiiten et al. 



20051 ). showing that environment is an important predictor of galaxy properties. 



Abundant in galaxy aggregations, clusters of galaxies offer the best sites to study 
such environmental effects. The observation that clusters at higher redshift contain a 
larger portion of blue galaxies is commonly referred to as the Butcher-Oemler effect (e.g.. 



Butcher fc Oemler 1984 



Poggianti et al.ll2006l : lEUingson et al.ll200ll : iGerke et al.ll2007n . The 



increase by ~ 20% in the blue galaxy fraction from 2=0. 03 to z=0.50 in the Butcher-Oemler 
effect suggests that star forming activity in galaxy clusters decreases toward the nearby 
Universe. With the continuous accretion of field galaxies into galaxy clusters, this decrease 
implies that star formation is inhibited by the cluster environment. 
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However, the fraction of galaxies located in these densest galaxy environments is only 
~ 10% even at the present epoch. Groups of galaxies, the intermediate density locations 
between dense clusters and the sparse field, are the most common environment where galax- 
ies are found. Seeded by lower-density perturbations in a hierarchical Universe, galaxy 
groups have a large variety of galaxy densities, from early-collapsed 'fossil groups' to aggre- 
gations of a few galaxies with densities only slightly greater than that of the field. With 
relatively small velocity dispersions compared with clusters, galaxy groups a re also the fa- 
vored environment for slow d ynamical encounters, leading to mergers (e.g., iMamonI l2000l : 
Coziol fc Plauchu-Fraynl 120071 ). Furthermore, a higher fraction of pa s sive galaxies have bee n 



observed in galaxy groups than in the field (e.g., IWilman et al.ll2005l . l2008l : lYan et al.ll2008l ). 
A sort of morphology segregation in galaxy groups is observed in the nearby Universe, as 
early-type galaxies are more commonly located in the g roup centers arid spi r als have larger 



proje cted centric distance from the host groups (e.g., IJeltema et al.l 120071 : IVerdugo et al. 



20081 ) . The underlying mechanisms which drive a sim ilar trend in the morphology-density re 



lation in clusters may also exist in galaxy groups (e.g.. lHelsdon fc Ponmanll2003al Jbl: lMulchaey fc Zabludofl 
1998h . 



A wide variety of observations in both optical and X-ray regimes has provided convincing 
evidence that substructures in clusters are comrnon. R i ch galaxy clusters a r e believed to be 



1996: Adami et al. 



built up by infalling gal a xy groups (e.g.. iBekkil Il999l : ISolanes et al.l Il999l : [Abraham et al. 



20051 : iKrick et all \20qS \. while galaxies in a group may have evolved 



already befo re the whole group falls into a cluster: i.e., what is referred to as pre-processing 
(lFujitall2004j ). However, there is still considerable uncertainty on the roles that these smaller 
structures play in the star formation history of cluster galaxies. Therefore, cluster and group 
galaxy populations give us a unique opportunity to study the environmental effect on large 
(global, or 'cluster') and small (local, or 'group') scales. 

In order to unfold the role of time and environment in galaxy evolution in clusters, 
we study these two major factors using both individual galaxies and galaxy groups in 16 
clusters from the Canadian Network for Observational Cosmology (CNOCl) multi-band 
BVRcI photometric observations. The structure of this paper is as follows. We present the 
data in |J21 We then briefly describe the photometric redshift method and sample selection in 
§31 The group finding method is briefly presented in §H The identification of cluster galaxies 
and groups is detailed in ^ We present the calculations of the environmental parameters 
and red galaxy fraction in §6] We describe the correction due to background contamination 
in f|71 We present the results ^ and discuss their implications in §21 Finally a summary is 
presented in §TU1 We adopt ^2^=0. 7, fiA=0.3, and Hq=70 km/s/Mpc. 
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2. The Cluster Catalogs 



2.1. The Data 



The original CNOCl Cluster Redshift Survey (lYee et al.lll996l ) contains 16 rich X-ray 
selected galaxy clusters at 0.17 < z < 0.55, and was carried out with the main goal of 
measuring the cosmological density parameter, VLm- The observations were made using the 
CFHT Multi- Object- Spectrograph (MOS) with Gunn g and r images to r ~ 24 and spec- 
troscopy to r ~ 21.5. The multi-color photometric follow-up observations in this study 
targeted 15 of the CNOCl clusters using the CFHT 12K camera in 2000 and 2001. The 
camera, composed of 12 chips of 2045 x 4096 pixels with 0.206" per pixel, has a 42' x 28' 
field of view. Ten clusters were taken from classical observations and the other five were from 
queue observations. Among these 15 clusters, 14 were observed using BVRcI filters with 
seeing ranging from ~ 0.70" to 1.24" in the i?c-passband. MS0015+16 (CL0015+16, z=0.55) 
was observed in B and V only, and hence is excluded for not having enough passbands for 
our photometric redshift method. We also exclude the MS1224+20 field because the obser- 
vations were non-photometric and not sufficiently deep for our scientific goals. MS1006+12 
was omitted due to a very bright star in the field. MS0302+16 {z=0A25, q;2ooo=03:05:35.42, 
52000=17:10:3.43) was targeted incorrectly to MS0302.5+1717 (;z=0.425, a2ooo=03:05:18.110, 
^2000= 17:28:24.92), about 18.82' north from the original CNOCl target. Another clus- 
ter, CL0303+1706 (z=0.418, a2ooo=03:06:18.7, (52ooo= +17:18:03), is in the south-east of 
MS0302.5+1717. 



1560 (a^non=12:34 :07.1, 6; 



2000- 



AbeU et al. 


1989; 


Ulmer & Cruddace 


1982 



David et al.lll999l ) is sou th of MS1231+15. M S1358-I-62 has been known as a binary cluster 
in the CNOCl analysis JCarlberg et al.lll996h . The other cluster in the MS1358+62 field is 
at 02000=13:59:39.511, ^2000=62:18:47.0, with 2;=0.329. In summary, our data set contains 
12 CNOCl clusters and 4 other confirmed clusters. The sample properties and exposure 
times are listed in Tables [1] and [2l 



2.2. Photometry and Color Calibration 



For the queue-observed data, the basic data reduction and calibration steps were un- 
dertaken using the Elixir system at CFHT 0. For the classically observed data, the data 



http://www.cfht.hawaii.edu/Instrunients/Elixir/realtinie. html 
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reduction proceeded in the standard way using IRAF (http://iraf.noao.edu/) including bias 
subtraction and flat fiel ding. Th e phot ometric calibration was performed using images of 
standard star fields from iLandoltl (Il992l ). 



The object finding, aperture photometry, and star-galaxy classification were carried out 
using the Pictur e Processing Package (PPP), details of which are described in lYed (119911 ) and 
Yee et al.l (119961 ). Here we present a brief summary of the procedures. This algorithm detects 
objects by the locations of brightness peaks above a preset threshold level and satisfying a 
number of criteria (e.g., a minimum number of connected pixels, a sharpness test, etc., 
see Yee et al. 1996). In general, about 2,000 to 3,000 objects on average were detected 
in each chip. The photometry is achieved by constructing and analyzing the fiux growth 
curve of an object. The total magnitude is measured in an 'optimal aperture', and then 
corrected to a standard aperture of 8.5" diameter using the PSF profile of reference stars, 
if the optimal aperture of the object is smaller than the standard aperture. The color of 
each object is measured independently from the total magnitude, using a 'color aperture' 
of 3", or the optimal aperture, whichever is smaller. The total magnitude for the second 
filter is then computed from the color and the total magnitude of the first filter. The star- 
galaxy classification is accomplished by comparing the growth curve of each object to a 
local set of bright, but not saturated, reference stars. The algorithm classifies objects into 
four categories: galaxies, stars, saturated stars, and spurious 'non-object' (e.g., cosmic ray 
detections or cosmetic defects). 

In order to check the photometry of the CNOCl follow-up observations, we compare our 
Rc data with the Gunn r photometry obtained by the original CNOCl program, using the 
transformation r = Rc + 0.28 + 0.027 x [B — Rc) derived from stars in the M67 field. To check 
the color offsets in the other passbands, we use the calibrated Rc photometry as the reference. 
For B and V p hotometry, we compar e the star colors with those in the R ed-Sequence Cluster 
Survey (RCS; iGladders fc Yed l2005l ) CFHT patches (IHsieh et al.l 120051 ) using bright stars 
(18 < Rc < 22). Because the RCS uses BVRcz' instead of BVRcI passbands, we compare 
the distributions of I — Rc star colors with those in the MS1231+15, MS1455+22, and 
MS1512+36 fields. These three pointings have been processed using the Elixir system, for 
which the /-band photometry is consistent with each other after being cross-checked. The 
zero points in BVI are shifted if any color offsets are present. 



2.3. Astrometry 



Astrometry for each frame is calibrated using the IRAF package MSCRED with 
USN0-a2 catalogs. The MSCRED package has the ability to operate on multi-chip CCD 
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data simultaneously without the need to split each image frame into separate CCD images 
and process them individually. Therefore, each CNOCl 12K follow-up observation pointing 
is made into a multi-extension FITS with 12 extensions. An astrometry calibration file is 
generated by the msctpeak command based on the MS0302-I-16 field. After updating the 
image headers with the calibration file, the mscczero and msccmatch commands are used 
to apply the shifts and match the positions between the objects on the images and in the 
USN0-a2 catalogs. We also update the astrometry in the headers of the original CNOCl 
MOS images, so the object list from the CFH-12K data can be easily matched with those 
spectroscopically observed using MOS. 



3. Photometric Redshift and Sample Selection 

The photometric redshift technique allows us to estimate galaxy redshifts using pho- 
tometry information. Because the information obtained from photometric redshifts is less 
accurate, the use of the technique requires careful consideration in selecting the samples. In 
this section we present briefly the photometric redshift technique and the sample selection; 
details can be found in Li & Yee (2008). 



3.1. Redshift Estimation 



Galaxy redshifts are estimated using an empirical photometric redshift technique from Li 
& Yee (2008). This technique assumes that redshift is a polynomial function of galaxy colors 
and magnitudes. The coefficients of this polynomial can be estimated by fitting the galaxies 
in a training set, which is a catalog containing galaxy spectroscopic redshifts and magnitudes. 
Our training set contains 398 8 galaxies in BVRrIz' uy to 2; ~ 1.4 with Rc < 24., constructe d 
from the Hubble Deep Field Jciavalisco et al.ll2004[) and CN 0C2 catal ogs JVee et al.l[200ol ). 
The properties of the training set are described in lHsieh et al.l (120051 ) andlLi fc Yeel J2008h . To 
achieve better-constrained solutions, the galaxies with spectroscopic redshifts in the CNOCl 
sample (2022 galaxies in total) are added to the training set. The photometric redshift 
uncertainties and probability densities of individual galaxies, P{z), are obtained as part of 
our photometric redshift technique as well. 

To assess the accuracy of our photometric redshifts, we compare them with the spec- 
troscopic redshifts of the 2022 CNOCl galaxies. The differences between photometric and 
spectroscopic redshifts are plotted in Figure [T] as functions of galaxy magnitude and B — Rc 



color. The overall r.m.s., a{zphot — Zspec), is ~ 0.051, with a{zphot — z, 



spec J 



0.048 and 
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~ 0.108, for Rc < 21.5 and Rc > 21.5, respectively. Blue galaxies (defined hj B — Rc < 1.8) 
has cr{zphot — Zspec) ~ 0.070, while a{zphot — ^spec) is ~ 0.041 for red galaxies (defined by 
B-Rc> 1.8). 



3.2. Galaxy Sample Selection 



For our analysis, we construct a sample of galaxies with a photometric redshift range 
of 0.02 < z < 1.4. The upper photometric redshift limit is due to the passband wavelength 
coverage for the 4000 A break in our training set. We select galaxies by their photomet- 
ric redshift probability densities. The photometric redshift probability densities used are 



consistent with the observed scatter in Zph, 



ot 



^spec 



in the training set for different types of 



galaxies. The selection criterion for galaxy i is set as: 



P,{z')dz' > 0.997, 



Zi—3<Jz 



with az set to 0.2(1+^;), which is equivalent to excluding galaxies whose photometric redshift 
uncertainty is greater than 0.2(l+z). 

Each galaxy in the sample is then assigned a completeness correction weight, Wi, to 
account for galaxies excluded from the sample due to the criterion. The completeness 
correction is computed using the ratio of the total number of galaxies to the number of 
galaxies satisfying our selection criteria in magnitude bins of Am/j^ = 0.1. We set a cutoff 
in Rc based on where Wi = 2 to avoid high weight galaxies. This apparent magnitude cutoff 
is brighter than the 100% photometric completeness magnitude in all cases, and ranges 
from 22.12 to 22.80 for the different pointings. The median and mean of the completeness 
correction factor Wi in the sample are ~ 1.08 and ~ 1.20, respectively. 

The sample used for group finding is limited to galaxies with k- and evolution-corrected 
Rr -hand absolute magnitude Mjf , brighter than M^^ + 2. We adopt M^^ =-21.41 at z = 
(IKodama fc Arimotd 119971 ). The k-correction is computed using the tabulated values 
from iPoggiantil (119971 ) for different spectral energy distributions (SED). For each galaxy, we 
interpolate the k-correction values between the E/SO and Sc SEDs based on its B — Rc color. 
We note that when selecting galaxies at different redshifts using k-corrections derived from 
galaxy colors creates a less biased sample th an applying k-c orrections averaged over SED 
types without using color information (e.g., see lLoh et al.ll2008l ). The Rc luminosity evolution 
is approx imated as Mf^:) = M(0) — zQ, where Q=1.24 for red galaxies and Q=0.11 for blue 
galaxies (ILin et al.lll999l ). In the cases when the apparent magnitude limit is brighter than 
M^^ -|- 2 at the photometric redshift of the group, we apply a correction factor R^ which 
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is the ratio of the integral of the galaxy luminosity function (LF) integrated to M^^ + 2 to 
that integrated to the completeness limit (see Eqn 2 in Li & Yee 2008), to correct for the 
expected missing galaxy counts in these fields. Note that while we use M^^ + 2.0 as the limit 
of pFoF group finding (which optimizes group-finding results), we limit the galaxy sample 
to M^^ + 1.5 for our galaxy population analysis. 



4. Identifying Galaxy Groups 

Galaxy groups in the photometric redshift catalogs are identified using the 'Probability 
Friends-of-Friends' algorithm (pFoF); a modified friends-of- friends group- finding algorithm 
which takes into consideration the photometric redshift probability density of each galaxy 
and grou p . Det ailed descriptions of the algorithm and performance tests are presented in 



Li fc Yed (120081 ): we provide a brief summary below. The main idea behind this algorithm 
is to treat the common photometric redshift space of the transversely connected galaxies as 
the group redshift space, which changes dynamically as new members are linked in. This 
updated group redshift is then used to search for more new members. Generally speaking, 
the group redshift is well confined by the group members, especially when the number of 
its members is large. The linking criteria are set by two parameters: DO^y and Pratio,crit- 
The DOxy is the reference physical transverse linking length at z=0, which is scaled by 
(1 + z)~^. In practice, the linking length includes the adjustments for the completeness 
correction, Wi, and the varying galaxy numbers at different redshifts due to the apparent 
magnitude cutoff. PraUo,crit is the criterion used in determining friendship in redshift space, 
where Pratio denotes the normalized total photometric redshift probability for a galaxy being 
in the same redshift space as another galaxy or a group of galaxies (see Li & Yee 2008 for 
details). The parameter Pratio essentially measures the amount of overlap in the photometric 
redshift probability distributions of the two objects being considered. Galaxies must have 
Pratio > Pratio,crit in Order to join into the group. We note that for the purpose of creating 
the group catalog the algorithm assigns a group membership to every galaxy , with isolated 
galaxies being identified as groups of one. 

This algorithm has been te sted using mock cat alogs constructed from the Virgo Con- 



sortium Millennium Simulation (jSpringel et al.ll2005l ). including the effects of using different 



linking criteria. With a fixed Pratio,crit, the recovery rate increases when larger DOxy is used; 
but the use of a larger DO^y also increases the fraction of falsely detected groups. With a 
given DOxy, the recovery rate is better when a smaller Pratio,crit is used, but the false detec- 
tion rate is larger as well. As a compromise between the recovery and false detection rates, 
the set of Pratio,crit = 0.37 and DO^y = 0.25 Mpc gives the best results. With these parame- 



-9- 



ters, the recovery rate is greater than 80% for mock groups of halo mass Mhaio greater than 
~ 1.2 X IO^^Mq, and the false detection rate is less than 10% when a pFoF group contains 
~ 8 or more net members (corresponding to M^aio ~ 3.7 x lO^^M©). The estimated pFoF 
group redshift uncertainty is ^Zgrp ~ 0.020. We adopt Pratio,crit — 0.37 and DO^y = 0.25 
Mpc for the group-finding procedure in this paper. 

For our analysis, wc select groups with richness Ngai > 8 and Ng^ > 7, where Ngai is 
the net weighted galaxy number, and Ng^ is the actual number of the linked group members 
returned by the pFoF algorithm. From the 13 pointings we identify a total of 188 groups 
satisfying our richness selection using the adopted pFoF parameters, including the cluster 
cores themselves. 



5. The Sample of Cluster Galaxies and Groups 

5.1. The Cluster as the Main Group 

Generally speaking, galaxy clusters can be treated as exceptionally rich galaxy groups. 
Because our sample is selected from the follow-up observations of CNOCl clusters with 
cluster redshift information available, we identify the galaxy group that coincides with the 
cluster center position as the main body, or core, of each cluster. 

The center of a group is determined using its members in the dense region where their 
local galaxy densities (see §6.1) are above the mean value of the group members. The center 
is calculated as the mean R.A. and Dec. of these members, weighted by their luminosity, 
completeness correction Wi, and local galaxy density. If a group's members are distributed 
in more than one dense region (i.e., clumps in the contours of local galaxy density), only 
galaxies in the clump which has the largest galaxy number counts are used. To designate 
the main galaxy aggregation, or the core, of a cluster, each galaxy group is assigned a score 
based on: (1) the inverse of the separation between the group and cluster centers, (2) the 
agreement of the group redshift compared with the cluster spectroscopic redshift, and (3) 
the group richness. The galaxy group which has the best score is chosen as the cluster main- 
group, i.e., the core galaxy aggregation of a cluster. We emphasize that cluster main-group 
associated with each galaxy cluster does not contain the whole of the cluster, but rather just 
the galaxies in the core that satisfy the hnking criteria of the pFoF algorithm. 

The redshifts of our 16 main cluster galaxy aggregations {zpFoF,ci) compared with the 
spectroscopic redshifts {zspec,ci) have a dispersion of cr{zpPoF,ci ~ Zspec,ci) ~ 0.019. The centers 
of our cluster main-groups are on average ~ 30" from the physical centers (i.e., where the 
cD galaxies are). The basic properties of our 16 identified cluster main-groups are listed in 
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Table O 



5.2. Cluster Galaxies and Groups 

To identify galaxies and groups in the cluster redshift space, we follow the same idea 
used in determining group membership in the pFoF algorithm (Li & Yee 2008), but applied 
only in the redshift direction without considering the transverse friends-of-friends criterion. 
That is, the PraUo of each galaxy or group with respect to the redshift probability density 
of the cluster main-group should satisfy Pratio > PraUo,crit, with the same PraUo,crit value 
used in the pFoF algorithm. The galaxies and groups that satisfy this condition are then 
considered to be in the large scale structure that make up the cluster, including galaxies at 
large cluster-centric radii. 

Applying the pFoF criteria to our group sample with respect to the cluster main-groups, 
we find 89 groups (not counting the cluster main-groups) consistent with being in the same 
redshift space as the clusters. Of these, 59 are within 3i?2oo of the cluster centers. Similarly, 
we apply the pFoF criteria to individual galaxies brighter than M^^ + 1.5 to create a sample 
of galaxies considered to be in the redshift space of the clusters. Hereafter, by 'cluster galax- 
ies' or 'cluster groups', we mean galaxies or groups whose photometric redshift probability 
densities within the cluster redshift space satisfy the Pratio,crit criterion. Based on the overall 
group sample and the typical uncertainty in the photometric redshifts for the groups, we 
estimate that the cluster group sample may be contaminated by projection (i.e., groups not 
associated with the same large scale structure in which the the galaxy cluster resides) at 
a rate of ~ 0.25 groups per cluster, sufficiently small that the effect can be ignored in our 
analysis. 

We note that two pointings, 0302-1-16 and MS0451-03, are not complete to M^^ + 1.5 at 
the cluster redshifts. For these clusters, we apply a correction factor, (see §3.2), based on 
the galaxy LF, to correct the counts to the M^^ + 1.5 limit. Note that this correction factor 
is not strictly correct when used for estimating the red galaxy fraction f^e A (see ji6.3), since 



the r ed and blue cluster galaxy LFs are not expected to be identical (e.g., iBarkhouse et al. 



20071 ). and will tend to decrease the estimated fred due to the steeper faint end of the blue 
galaxy LF. However, since our adopted absolute magnitude limit is relatively bright and 
the incompleteness ranges from 0.05 to 0.3 mag, this effect does not significantly affect our 
conclusions. 



We present an example in Figure [2] using the cluster Abell 2390 {z = 0.228) to illustrate 
the selection of cluster galaxies using photometric redshifts. The cluster galaxies, whose 
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memberships are determined using the redshift assigned by the group finding algorithm, 
exhibit an excess peak in the photometric redshift space in contrast to the field galaxy 
distribution. Our method of selecting cluster galaxies and groups with reassignment of 
their redshifts to those of the associated pFoF groups eliminates foreground and background 
galaxies with higher efficiency than the usual method of applying a photometric redshift cut 
based on the individual galaxy's original photometric redshift and uncertainty. 

Figure [3] shows the sky locations for the cluster galaxies and cluster groups in the Abell 
2390 field. The cluster has an elongation fro m the central core at a bout -50° from th e N-S 
axis, spanning both sides of the cD galaxy (lAbraham et al.l Il996l : iPierre et al.lll996l ). An 
infalling s ub-component at the projected distance of ~ 650" to the cluster center has been 
studied in lAbraham et al.l (Il996l ). This group is identified by the pFoF algorithm as well, and 
is marked by crosses in Figure [31 Our pFoF algorithm also identifies two additional groups 
in the same direction, which are marked as squares and triangles in Figure [31 However, as 
they lie about 100" and 280" away from the spectroscopically identified group, these two 
groups do not contain any members in the original CNOCl spectroscopic coverage. 

We also present the observed color-magnitude diagrams in Figure H] for all these cluster 
groups. Galaxies in the cluster main-group exhibit a well-defined red sequence, while galaxies 
in the other groups have visually somewhat less prominent red sequences, primarily due to 
the small number of galaxies. We note that galaxies in groups in clusters may also include 
some galaxies from the main cluster in projection, which will produce systematically redder 
groups closer to the cluster core. However, this effect on the average colors of the richer 
cluster groups used in our analysis (§8.4) is likely to be minimal. The typical galaxy surface 
density for groups with 8 members or more is 21.56 galaxies/ Mpc^, whereas the average non- 
group cluster galaxy surface density at r(7i=0.5-1.5 ranges from 2.43 to 5.95 galaxies/Mpc^, 
and is as low as 2.05 galaxies/Mpc^ at the largest radii we probe. 



6. The Environmental Parameters and Red Fraction, fred 

We intend to study the colors of cluster galaxies in different environments with the aim 
of identifying environmental effects on galaxy evolution. For this purpose, we demarcate 
the environmental effect into local and global regimes, and use the fraction of red galaxies as 
an indicator of their evolutionary stage. Here, we present our methods for computing the 
environmental parameters and the red galaxy fraction. 
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6.1. Local Galaxy Environment: Local Galaxy Density E5 

For each galaxy, the local galaxy density, S5, is computed using a circular aperture and 
counting galaxies to < M^^ + 1.5: 

= ^ ^bgizci)- (2) 

In this calculation, n is the n^^ nearest galaxy at a distance r„ from the seed galaxy such 
that the total completeness weight Wi > 5. When there is a non-unity weight correction, 
rs is not necessarily the distance r„, since the total number of objects at r„ would not be the 
integer 5. Assuming a uniform density, we can approximately correct for the effect of the 
'missing galaxies' represented by the weight corrections by computing the in Equation [2] as 

^5 = \l^^^^^n- The computation of S5 includes a background correction term T,iig{zcL) 
(see §7). In Figure [3l we also overlay the contours of local galaxy density to illustrate the 
substructures in the cluster redshift space for Abell 2390. 



6.2. Global Cluster Environment: Cluster-Centric Radius tcl 



The global cluster environmental influence is delineated using the cluster-centric radius 
parameter, denoted as rcL, in units of -R200, where i?2oo is the radius within which the mass 
density is 200 times the critical density. We use the i?2oo values published in lCarlberg et al. 



(I1997I ) for our CNOCl clusters, except for MS0906+11. Note that these are based on velocity 
dispersion measurements. The i?2oo value for MS0906+11 is overestimated, due to the binary 
nature of the cluste r, with two velocity structures appearing to be overlaid on the sky. 
Borgani et al.l (119991 ) also analyze the velocity dispersions of CNOCl clusters. In their work, 
they separate the two merging clusters in MS0906+11 and present their respectiv ely velocity 
dispersions as 886 and 725 km/s. Scaling with the value in ICarlberg et al.l (119971 ). we adopt 
i?2oo=l-79 Mpc for MS0906+11. For those clusters with unknown -R20O) we estimate the 
values using a solution derived from the correlation between Ngai and -R200: Ngai ~ 11-8 x 
i^^QQ. The values for the clusters are listed in Table [3l 



6.3. Red Galaxy Fraction: fred 

In a given local-density environment or within a galaxy group, the red galaxy fraction, 
fred, is computed using galaxies with M^'^ < M^^ + 1.5. We compute fred using the observed 
B — Rc color- magnitude diagram. First, an offset in the theoretical red-sequence zeropoint 
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(from iKodama fc Arimotd (119971 )) is applied based on our observed cluster data. We then 
take the color half way between the SEDs of E/SO and Sc as the boundary separating the 
red and blue galaxies. This boundary is close to the valley in the color distribution between 
the red sequence and the blue cloud. It ranges from ~0.39 mag at z = 0.2 to 0.63 mag at 
z = 0.5 in B — Rc to the blue of the red sequence. 

Background galaxies play an important role in the estimation of the true fred, especially 
in environments where the net galaxy counts are small. We therefore use Bayesia n statistical 
inference based on Poisson statistics to estimate the fred probability function (ID'Agostini 
2004J : lAndreon et al.ll2006l ). We then take the mean and the 68% confidence interval of a 
fred probability function as the estimated fred value and its uncertainty. 



7. The Control Sample 

When photometric redshifts are used to identify cluster galaxies, the background galaxies 
still have a significant role in contaminating the sample. Such contamination, however, can 
be statistically estimated and corrected fo r using a large contro l sample. We us e three 



RCS CFHT patches {0920, 1417, and 1614; iGladders fc Yeell2005l : iHsieh et al.ll2005l ) as the 



control fields. These three patches have the deepest photometry in the RCS survey and 
cover a total of 11.83 square degrees. To derive the galaxy surface density as a function of 
redshift, we first count the number of galaxies at each redshift, N{z), within a bin size of 
0.01 in z by summing all the weighted photometric redshift probability densities within the 
redshift interval: 

/•2+0.005 

N{z)= / P^{z')w4z', (3) 

^2-0.005 

where Pi{z) is the photometric redshift probability density of each galaxy and Wi is the 
complete weight (see §3.2). The surface density at each redshift n{z) is accordingly obtained 
by dividing N{z) by the total area of the RCS control samples at each redshift. The T,bg{zcL) 
in Eq. [2] is therefore computed as: 

POO 

^bgizci) = / L{zcL)n{zcL)dzcL, (4) 
Jo 

where L{zcl) is the photometric redshift likelihood for which the members of a cluster 
main-group may occur (see Li & Yee 2008). 

When we conduct analyses using subsamples of cluster galaxies based on the local galaxy 
density, the background corrections are also derived using galaxies selected from the control 
sample with the same criteria. For example, to estimate Nbg within a E5 region, the n{z) 
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in Eq. H] is computed using background galaxies with S5 which are smaller than the upper 
boundary of the S5 bin, instead of using only background galaxies whose S5 are within the 
same S5 range. This is because cluster galaxies in a high S5 region may be contaminated 
by background galaxies from both high and low-density regions, and those in low E5 regions 
are only affected by background galaxies in low E5 bins (otherwise their S5 would be larger). 
The area of a E5 clump, whe re S5 is above a fixed cutoff, is computed using the same method 
in calculating group area in iLi fc Yed (120081 ) . 



8. Results 



8.1. The Butcher-Oemler Effect 



The Butcher-Oemler effect (iButcher &: Oemlerl Il984l ) states that galaxy clusters at 
higher redshifts possess a larger fraction of blue members than those at lower redshifts. 
Instead of using the fraction of blue galaxies as the parameter to indicate the evolutionary 
stage of galaxies, we adopt the red galaxy fraction fred, where a larger fred reflects a smaller 
blue galaxy fraction and vice versa. T he Butcher-Oeraler effect is suggested to have a de 
pend ence on cluster centric radius (e.g., lEUingson et al.ll200ll : iDahlen et al.ll2004 iLoh et al. 
20081 ) ■ so we compute the fred for each cluster by using (1) members of each cluster main- 
group (as deflned in §5.1), (2) all cluster galaxies within I-R2005 and (3), within l-1.5i?2oo- 
We note that the cluster main-group members on average occupy an area smaller by a factor 
~ 4 than that within I-R200; i-e., an area with an equivalent radius of ~ 0.5i?2oo- The fred 
are computed using galaxies brighter than M^^ + 1.5. We note that this limit is identical to 
the original luminosity limit of My = —20 used by Butcher & Oemler (1984) after correcting 
for the different Hq used. 

In Figure Owe plot the fred of each cluster versus its spectroscopic redshift for different 
cluster-centric radii. We also overlay the results of linear fltting between fred and redshift. 
Figure [5] shows that there is a dependence of fred on the cluster-centric radius. The fred 
computed using galaxies in the cluster main-group have the largest values, and this subsample 
also exhibits the most gentle gradient . The fred for all three counting radii show a decreasing 
trend with redshift. It was shown in lEUingson et al.l (1200 ll ) that the Butcher-Oemler effect 
is not signiflcant within 0.5i?2oo; our result of Figure [5] is consistent with their work. For 
z < 0.40, fred is essentially constant for the cluster main-groups (i.e., in the cluster core), 
but at higher redshift there appears to be a drop in the distribution of fred, with two out 
of the four clusters at ^ > 0.40 having fred signiflcantly lower than the values from the low- 
redshift clusters. This result can be interpreted as the fred in the cores of clusters reaching 
a 'saturation' value close to 1 at z ~ 0.4, while the Butcher-Oemler effect continues to be 
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observed further out in radius. 



We note that the fred in the core of MS0451-03, at z = 0.539, is a significant outher, 
in that its value is comparable to those at larger radii. This is likely a stochastic effect 
due to cosmic variance and not a general redshift trend, as the other CNOCl cluster at the 
same redshift, MS0016+16 (for which we do not have photometric redshift data), shows a 
considerably higher red galaxy frac tion in the core, deriv ed using the CNOCl spectroscopic 
data, than that of MS0451-03 (see lEUingson et al.ll200ll ). A sudden change in the average 
core red fraction s of clusters is not seen in the large sample of clusters at 0.45 < z < 0.90 by 
Loh et al.l (120081 ) . further suggesting that MS0451-03 is significantly bluer in the core than 
the average for clusters at this redshift. 



8.2. The Butcher-Oemler Effect and Stellar-Mass Selected Samples 



The Butcher-Oemler effect has 
in 1984 (e.g., 



Margoniner et al 



20051 : iLoh et al.l l2008l : iNuiiten et al. 



j rought much debate i n the literature since its first study 



2001 



Metevier et al.ll2000l : IPe Propris et al.ll2004l : iTran et al. 



20051 ). The most controversial one is its existence: 



whether the Butcher-Oemler effect is a real phenomenon or simply due to sample selection 



effects (e.g., lAndreon &: Ettorilll999l : lAndreon et al.l 120041 ). For example, there are studies 



showing that the Butcher-Oemler effe ct becomes insignifican t when the gal axy sample is 



selecte d using K-band luminosity (e.g.. lDe Propris et al.ll2003l ). More recently, iHolden et al. 



( 120071 ) . using a small sample of 5 rich clusters with 0.02 < z < 0.9 and sampling with 
a relatively high stellar-mass limit of 10^^'^ Mq, concluded that while the Butcher-Oemler 
effect is seen in a luminosity-selected cluster galaxy sample, a stellar-mass selected sample 
show insignificant changes in the galaxy red fraction over the redshift range. 

To examine this issue, we investigate the effect of stellar-mass selection on the Butcher- 
Oemler effect for the CNOCl sample. We derive the stellar r nass in sola r units , , for galax- 



i es in the cluster galaxy sample by applying the algorithm in iBell et al.l (120031 ) . iHolden et al. 



(120071 ) used the same method, and provided a detailed discussion of the robustness and uncer- 
tainties in t he derived stellar mass. We use the relation: log(M^,/L/?^) = -0.523 -|- 0.6 86(Mb 
- Mpic) from lBell et al.l (120031 ) with the values of the solar units from IWortheyl (119941 ). Both 
Mjic and Mb — M^c in our computation are k-corrected using color- dependent k-corrections 
(see §3.21) . We note that u s ing th e redder i?c-band luminosity, compared to the B band which 
was used by IHolden et al.l (120071 ). produces more robust mass estimates with less scatter. As 



examples, we present in Figure [6] the relationship between M'^^ and stellar mass for Abell 



2390 and MS1358-I-62. In the Figure we also show lines marking the sampling limits 
of M^^ + 1.0 and M^^ + 1.5, and stellar mass cuts at log = 10.6, 10.2, and 9.75. We 
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note that the mass cuts at lO^'''^ and 10^'^^ Mq correspond approximately to the average 
stellar mass for red and blue galaxies at ~ M|j^ + 1.5, respectively; whereas a stellar mass of 
IO^^'^Mq coincides with red galaxies at ~ M^^ + 1.0. Other clusters have remarkably sim- 
ilar plots, with the major difference being the completeness level (with most lower redshift 
clusters being considerably more complete than M^^ + 1.5). The most incomplete cluster is 
MS0451-03, which is complete to about M^^ + 1.2. We note that our Wi and Rw correction 
factors partially alleviate this problem. 

Using the stellar- mass selected sample for galaxies within I-R2005 we re-derive fred using a 
statistical background subtraction method similar to that for the luminosity-selected sample. 
Note that because of the considerably larger uncertainties in the photometric redshifts of 
the control sample galaxies, the derived stellar-mass measurements are expected to have 
significantly larger scatter. We show the results in Figure [TJ where we plot the red galaxy 
fractions for the three different stellar-mass cuts. All three samples show the Butcher- 
Oemler effect, with the lower mass limit sample showing the largest change in f^ed with 
redshift. We note that there is significant incompleteness in the red galaxies for the higher 
redshift clusters for the lowest mass cut (logM^, = 9.75) where our incompleteness correction 
procedure may not be effective and hence may produce smaller f^ed for the higher-redshift 
clusters. However, the log M^, = 10.2 sample, which is largely complete, also shows that a 
lower mass cut produces a more significant Butcher-Oemler effect. The result with the log 
M^, > 10.2 sample is essentially identical to that in Figure [5] (open diamonds, using galaxies 
inside I-R200 with a luminosity limit of ~ M|j^ + 1.5), except for a systematic shift to slightly 
higher average f^ed values, which is due to the larger number of red galaxies in the sample 
when a stellar-mass cut is used. 

Our sample shows that the slope of the fred-z trend is significantly shallower when using 
a higher stellar-mass limit. This is not a surprising conclusion, and is expected from the 
well-established observations pointing t o the down-sizing scenario of galaxy evolution (e.g.. 



up of the faint end of the red sequence with decreasing redshift (e.g., 


Gilbank et a 


2008 


van den Bosch et al. 


2008; 


De Lucia et al. 


2007; 


Tanaka et al. 


2005: 


De Lucia et al. 


2004) 



Both of these phenomena would naturally reduce the observed Butcher-Oemler effect, if the 
chosen sample has too high a limit in stellar mass or luminosity. Hence, the lack of strong 
evolution in the iHolden et al.l (120071 ) result may arise from their relatively high stellar-mass 
limit. The small sample size of 5 clusters used by Holden et al. also makes detecting a reduced 
Butcher-Oemler effect more difficult. Furthermore, the Holden et al. results are based on 
cluster galaxies within a cluster-centric radius of 1.25Mpc, likely significantly smaller than 
the i?2oo radii of these rich clusters. This will also further reduce the detectability of the 
Butcher-Oemler effect. 
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Figure[6]also demonstrates that a sample chosen using a SED-dependent k- and evolution- 
corrected luminosity limit can in fact provide an unbiased measurement of the Butcher- 
Oemler effect, if this effect is generalized to mean the evolution of the color fraction of 
galaxies. For a chosen M'^^ limit, there would be different limits for the average stellar mass 
for the red sequence galaxies and for the blue galaxies. As long as these mass limits do not 
change as a function of redshift, an M^^*^-selected sample essentially samples similar galaxies 
over the redshift range, and hence the evolution of fred is not driven by selection effects. 
As an example, for the CNOCl sample presented here, our k- and evolution-corrected M'^^ 
limit samples red and blue galaxies to stellar-mass limits of log of ~ 10.2 and 9.75, re- 
spectively, for all our clusters. Figure [8] plots the median stellar mass in the magnitude bin 

+ 0-5 < M'^I^ < M^^ + 1.5 for the red and blue cluster galaxy samples as a function 
of cluster redshift. This illustrates that there is no substantial selection effect in the stellar 
mass sampled as a function of redshift. The somewhat higher mass value for MS0451-03 is 
due to incompleteness in M^^ in the highest redshift cluster. 

In the remainder of the paper, we analyze the cluster galaxy stellar population using the 
luminosity-selected sample defined in §3.21 which is equivalent to sampling mass limits of log 
M,, = 10.2, and 9.75 for red and blue galaxies, respectively. This allows for a more robust 
background subtraction based on observed magnitudes, and also increases the number of 
galaxies used in deriving the various statistics. 



8.3. Cluster Galaxies 



In this subsection, we examine the dependence of galaxy populations on redshift, vcl, 
and S5 using samples of galaxies combined from the different clusters. The cluster galaxies 
are separated into three redshift bins: 0.15 < 2; < 0.30, 0.30 < z < 0.40, and 0.40 < z < 0.55, 
with 8, 4, and 4 clusters in each bin, respectively. 

We plot in Figure M the radial dependence of fred at each redshift bin using galaxies in 
all local density environments. As expected, there is a strong radial dependence of fred on 
cluster-centric radius for all redshift bins. Furthermore, we see a Butcher-Oemler effect at all 
radii between the high- and low-redshift bins. Figure [TO] shows fred as a function of E5 in each 
redshift bin using all cluster galaxies within rcL < 3-R2oo- The galaxy red fraction in general 
increases with E5. This is consistent with the well-established observational trend that that 
early-type, non-star-forming galaxies tend to populate hi gh S5 regions, while late-type, star- 



forming galaxies ar e more common in low S5 regions (e.g..lDresslerlll980l:lKodama et al.ll2001 



Cooper et al.l 120071 : iHolden et al.l 120071 : iGomez et al.ll2003l : iBlanton et al.ll2005h . Figure [TO] 
also shows that the magnitude of the Butcher-Oemler effect depends on the local galaxy 
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density. The high-density regions show a stronger increase in fred as a function of redshift, 
with the highest E5 bins reaching the nominal 'saturation' fred values close to 1 at z ~ 0.40; 
whereas the lower density bins show a progressive increase in fred over the whole redshift 
range. The trend in Figure [10] can be interpreted as galaxies in higher and higher density 
regions completing their star formation (or have their star formation quenched) at earlier 
and earlier epochs. For example, galaxies with S5 > 100 reach fred ~ 0.9 by z ~ 0.35, while 
those with 50 < S5 < 100 become mostly red the lower redshift of ~ 0.20. 

Because the apparent effects of tcl and S5 are largely degenerate, we probe the change 
in fred in different environments in more detail by computing its values in bins of S5 within 
a fixed rci location, and as a function of redshift. We divide the data into bins with 
rcL < 0.50, 0.50 < r^L < 1-0, 1.0 < r^L < 1-5, and 1.5 < r^L < 3.0. The data are also 
separated into bins of < S5 < 10, 10 < S5 < 20, 20 < S5 < 40, and E5 > 40. Figure [H] 
presents the 'fred-z' trend in each rci and S5 bin. In all but one panel, fred declines with 
increasing redshift, indicating that within the uncertainties of the data the Butcher- Oemler 
effect is observed over the whole ranges of tcl and S5 values. At large cluster-centric radii, 
we see that the Butcher-Oemler effect becomes progressively steeper with increasing local 
galaxy density. In contrast, in the cluster core the galaxy populations are dominated by red 
galaxies over most redshift bins, and show little evolution for z < 0.35 for all local densities. 

We examine the effect of local galaxy density on fred further in Figure [121 where we plot 
the fred as a function of E5 (the 'fred-^5 ' trend) in different tcl and redshift bins. In general 
we find an increase of fred with increasing S5, as expected from previous studies. However, 
the magnitude of the effect appears to be weak, and dependent on both the cluster-centric 
radius and the redshift. At each redshift, the slope of the '/^ed-Ss' trend is steepest at the 
outskirts (1.5 < tcl < 3.0), and very weak or non-existent for galaxies within the virial 
radius of the clusters. An interesting observation is that the overall increase in fred, as one 
moves from the outer radial bins to the inner ones, is greater than, or comparable to, the 
increase in fred as a function of ZI5 over the entire local galaxy density range. This, along 
with the general lack of a strong dependence of fred on S5 in the inner parts of the clusters, 
implies that the influence of the global cluster environmental, rather than the local galaxy 
density effect, dominates inside the cluster virialized region. Compared with the strong 'fred- 
S5 ' trends shown in Figure [TUl the weak correlations of fred with S5 seen when the data 
are divided into rci bins (Figure [12]) suggest that at least part of the fred-'^5 dependence in 
Figure [TO] is the result of a cluster-centric radius effect, since galaxies in high-density regions 
have a higher chance of being found in small tcl locations. 

In Figure [13], we plot fred as a function of rcL in different redshift bins to examine 
more explicitly the effect of the global cluster environment on galaxies, with the local galaxy 
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density controlled. We find that, in general, fred decreases with increasing tcl for all local 
densities. For the two lower-redshift samples, high S5 regions have, on average, flatter 'fred- 
TcL ' trends than those for the low local density bins. However, the 'fred-^cL ' trends for the 
high-redshift sample are relatively steep for all the density bins. This indicates that, over 
our redshift range, the apparent effect of the global cluster environment is strongest in low 
local-density regions. This result suggests that galaxies in high local-density regions (e.g., 
galaxy groups) are likely to be already in a more advanced evolutionary stage, i.e., with large 
fred, before infalling into the cluster, and hence the global effect of the cluster on their fred 
would be less apparent. 



8.4. Cluster Galaxies in Galaxy Groups 

From the results of §8.3[ it is apparent that cluster galaxies in high local density re- 
gions have a different evolutionary history from those in low-density environments; however, 
there is also evidence that the observed galaxy density may not be the primary parameter 
producing these effects. High galaxy density regions outside the cluster core are often indi- 
cations of galaxy groups; therefore, it is useful to attempt to delineate the effects of a group 
environment from those of a high local galaxy density environment. 

To probe the role of galaxy groups in the evolution of cluster galaxies, we flag cluster 
galaxies into 'group' and 'non-group' cluster galaxies according to the pFoF group flnding 
results (see § 15. 2p . Group cluster galaxies are the cluster galaxies classifled as the members 
of pFoF groups of Ngai > 8 and Ng^ > 7, and non-group cluster galaxies are those considered 
by the group- flnding algorithm to be in pFoF groups with Ngai < 2 and Ng^ < 2. Recall that 
the pFoF algorithm assigns a "group" membership to all galaxies, including isolated galaxies 
(which are assigned to groups with one member). The gaps in the Ngai and Ngz cutoffs serve 
to reduce the ambiguity in the group membership of the galaxies. The relatively high Ngai 
criterion, which corresponds to a halo mass of ~ 3.7 x IO^^Mq, is chosen so that the group 
sample has an acceptable false detection rate of < 10%. The low Ngai criterion ensures that 
the non-group cluster galaxies are relatively isolated and in subhalos with no more than two 
galaxies brighter than M^^ + 2.0. Scaling with Ngai, we expected these galaxies to be in 
halos with mass less than lO^^M©. The pFoF groups with 3 to 7 members (the 'in-between' 
sample) are less well deflned, as they have high false detection rates, ranging from 70% to 
20%. These groups are expected to have a broad halo mass distribution, from that of a single 
galaxy to ~ 3 x lO^^M©. The cluster galaxies associated with the cluster main-groups are 
excluded from the group cluster galaxies, so that the environmental effects caused by galaxy 
groups can be speciflcally investigated. We can consider that cluster main-group galaxies 
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are the galaxies residing in the core of the main dark matter halo which defines the cluster. 

We have a total of 1073, 1226 and 3507 cluster galaxies in the cluster main-group, group, 
and non-group subsamples, respectively. A total of 5727 galaxies fall in between the group 
and non-group criteria. It is interesting to note that with this subdivision of cluster galaxies, 
we find that ~ 17.5% of the galaxies outside the virial radius (> I-R200) are in subhalos 
with masses greate r ~ 3-4 x IO^^Mq. This is broadly consistent with the simulations of 



Berrier et al.l (120091 ) which found ~15% of galaxies in their most massive clusters (which are 
most comparable with the CNOCl sample) are accreted from groups with halo masses larger 
than ~ 3.5 x IO^^Mq. 

We illustrate in Figure [H] an example of the distributions of the three subsamples of 
cluster galaxies in the J^^-tcl plane using Abell 2390. For completeness, the 'in-between' 
sample (between those considered to be group cluster and non-group cluster galaxies) is also 
plotted. The S5 distribution can be described in general as having a decreasing envelope 
with increasing radius, with high-density groups super- impose d on it. The envelop e distri- 



bution can be fitted reasonably well with a projected NFW (iNavarro et al.l Il996l ) profile. 
It is interesting to note that Figure [H] demonstrates that the E5 measurements computed 
using photometric redshifts can produce a reasonable projected galaxy density profile for 
the cluster, indicating that these density measurements are in fact physically meaningful, at 
least when considered over a factor of several in density. 

Figure [13 shows the redshift dependence of fred for the three subsamples of galaxies 
inside rcL=3, using galaxies in all local galaxy density environments. We see that all three 
samples show an increase of fred with decreasing redshift, but have systematically different 
levels of fred, with the cluster main-group being the reddest and the non-group cluster 
galaxies the bluest. Furthermore, the three subsamples also show different evolution in 
fred- The cluster main-group galaxies exhibit the largest fred and reach the final state of 
fred ~ 1 earlier than the group cluster galaxies. The group cluster galaxies show the strongest 
evolution from z ~ 0.5 to 0.2, reaching fred close to that of the cluster main-group galaxies 
at the low-redshift bin. The non-group cluster galaxies, located in less massive subhalos, 
show a milder increase in fred, still retaining a substantial fraction of blue galaxies in the 
low-redshift bin. They begin to partake in significant evolution at the lower redshift of 
z ~ 0.35. In the context of the halo model, the cluster cores are the center of the most 
massive dark matter halos, while group and non-group cluster galaxies are in successively 
less massive subhalos. The epoch when fred reaches ~ 1 is therefore a refiection of the 
dependence of galaxy evolution on the halo mass: galaxies in more massive halos have their 
star formation truncated at an earlier time than do those in less massive halos. This result 
can be described as a 'group down-sizing' effect, akin to the down-sizing effect seen in the 
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evolu tion of individual galaxies (e.g.. lCowie et al.lll996l : ICattaneo et al.ll2008l : ISeymour et al. 



20081 ). though the physical mechanisms underlying this effect may be very different. 



Figure [12] demonstrates the different evolutionary histories of galaxies in subhalos of 
different sizes by plotting the dependence of fred on cluster centric radius for the three 
galaxy subsamples in different redshift bins. We plot the cluster main-groups as a single 
point at a small radius. Cluster galaxies in groups show the largest increases in fred with 
decreasing redshift; and this increase is seen in both inner and outer cluster-centric radius 
bins. The fred for galaxies in groups in the cluster outskirts shows the largest change with 
redshift, catching up to the value of the cluster main-group galaxies in the core by z ~ 0.2, 
by which time a fiat fred distribution in cluster-centric radius is seen. In comparison, the 
fred of non-group cluster galaxies shows a very strong radial gradient at the low-redshift bin. 
The rate of evolution of fred for the non-group cluster galaxies depends significantly on the 
cluster-centric radius, such that galaxies in the outermost radial bin (at ~ 2i?2oo) show little 
change with redshift, while galaxies near i?2oo evolve significantly. 

Our results suggest that galaxies in moderately massive groups outside of the cluster 
core are turning red in response to the environment of their own dark matter subhalos, with 
the cluster environment having a strong effect possibly only at the core. On the other hand, 
the non-group cluster galaxies, situated in low mass subhalos, appear to turn red at a low or 
non-existent rate at the outskirts, and significant evolution is seen only when they are near 
the virial radius of the cluster. This result can be interpreted as galaxies in small subhalos 
being much more affected by the massive cluster dark matter halo into which they are falling. 
However, another possible contribution to this increase of fred for the non-group sample near 
the virial radius could come from member galaxies of more massive infalling groups (which 
have already turned red), which have been dynamically stripped from their parent subhalos. 

To separate the effects of local galaxy density and group membership, we plot in Figure 
[T71 fred for the three subsamples at different redshifts as a function of S5. The uncertain- 
ties are fairly large due to the small sample size, but, nevertheless, the plots suggest that, 
after controlling for the local galaxy density, the three subsamples still have different red 
fractions. This result indicates that the different average local galaxy densities in these 
subsamples do not entirely explain the different fred measured. Comparing galaxies in the 
cluster main-groups to those in groups in the high-redshift bin, the cluster galaxies in groups 
have substantially smaller fred at all densities, similar to non-group cluster galaxies. By 
z ~ 0.2, the group cluster galaxies have evolved into having fred similar to the cluster main- 
group galaxies, while the non-group cluster galaxies show a considerably milder evolution. 
Thus, by 2 ~ 0.2, cluster galaxies in groups are redder than their non-group counterparts in 
similar local galaxy densities, implying that group membership, or more precisely, the mass 
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of the subhalo in which a galaxy resides, has a stronger effect on the evolution of its stellar 
population than does the local galaxy density. 

9. Discussion 

The availability of photometric redshift information over a large field of view allows us 
to examine the Butcher-Oemler effect in galaxy clusters in detail, including investigating the 
effects of local galaxy density and the dependence on the global cluster scale. Furthermore, 
by separating the cluster galaxies into different group status categories, roughly indicative 
of the size of the dark matter subhalo in which they reside, we are able to examine more 
closely the nature of the dependence of the galaxy red fraction on their subhalo masses. 

Figures [9] to [13] illustrate the apparent dependence of galaxy colors on local galaxy 
density and cluster centric radius, averaged over the whole cluster galaxy sample. Looking 
at cluster galaxies as a whole, by z ~ 0.20, most cluster galaxies in high S5 regions already 
have their star formation truncated, as their fred is close to 1 (see Figure 10). The truncation 
process appears to occur earlier for galaxies in denser regions. Assuming a galaxy would turn 
red 1~2 Gyr after star formation is quenched, this truncation process is mostly completed 
by the redshift of ~ 0.45 for galaxies in high-density regions, given that their fred reaches 
~1 by our low-redshift bin {z ~ 0.20). 

The cluster environment, as defined by the position of the galaxy in the cluster, however, 
appears to have an even stronger effect on f^ed than the local galaxy density. For galaxies with 
similar E5, those closer to the cluster center have a larger fraction of their star formation 
suppressed than those in the outskirts (Figure 11, comparing within each row). When 
the galaxy sample is separated into subsamples in cluster-centric radius, the dependence of 
fred on the local galaxy density is considerably weaker, and nearly non-existent in the cores of 
clusters (Figure 12), suggesting that local galaxy density is not the dominant environmental 
parameter. We also find galaxies in low-density regions have a steeper fred-^cL gradient, 
compared to those in high-density regions (Figure 13). This suggests that the effects of 
the global cluster environment are stronger for galaxies in low-density regions. However, an 
alternative interpretation might be that star formation in galaxies in high S5 regions may 
have already been largely stopped by some other processes not associated with the global 
cluster environment; whereas in low S5 regions, there is still a relatively large fraction of star- 
forming galaxies available to have their star formation truncated by the cluster environment 
as they approach the cluster virial region. 

The weak dependence of fred on local galaxy density seen when the data are analyzed 
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in different tcl bins (Figure [T2|) suggests that S5 may not be a fundamental parameter in 
determining tlie evolution of the galaxy population in clusters. Separating cluster galaxies 
into group and non-group galaxies, representing galaxies in more and less massive subhalos, 
provides us with additional insights into cluster galaxy population evolution, and a more 
physical interpretation of the apparent dependence on local galaxy density. The dependence 
of fred on S5 for both group and non-group cluster galaxies is also weak (Figure [T71). similar 
to that found for the (whole) cluster galaxy sample after they are separated into bins of 
cluster-centric radii (Figure [T^ . The most significant result from applying the group and 
non-group selection of galaxy subsamples is presented in Figure [T6l It presents clear evidence 
that, outside the cluster core, galaxies situated in the more massive subhalos (the group 
cluster galaxies) have a much stronger Butcher-Oemler effect, regardless of their position 
relative to the cluster center. In the context of a dark matter halo model, this suggests that 
much of the dependence of the galaxy population on local galaxy density can be attributed 
to whether an inf ailing cluster galaxy is in a massive subhalo. 

Our analysis shows that cluster galaxies in groups, regardless of their local galaxy density 
and cluster-centric radius, exhibit a much stronger Butcher-Oemler effect than the non-group 
cluster galaxies (FigurefTBlandfTTI): the difference between these subsamples is especially large 
beyond the virial radius. In fact, the cluster galaxies considered to be in moderately massive 
subhalos appear to be the most significant contributors to the observed Butcher-Oemler 
effect at the cluster outskirts (~ 2 — 3-R200) over the redshift range of our sample. This 
dependence on halo mass can also be extended to include the core of a cluster, which is the 
center of the most massive halo. The cores of the clusters have high red fractions close to 
1 at 2; ~ 0.4; they complete the quenching of star formation the earliest, ahead of the less 
massive infalling subhalos. 

The evolution of f^ed of the large sample of cluster galaxies with pFoF group membership 
falling in between the group cluster and non-group cluster galaxies is also consistent with 
this picture. This sample contains the bulk of the cluster galaxies - about 49% of the cluster 
galaxies outside the cluster core. The fred of this sample at the outer cluster-centric radial 
bins show a Butcher-Oemler effect that is between those of the group cluster and non-group 
cluster galaxies. For example, at the outermost tcl bin, fred rises from 0.19 ±0.02 at z ~ 0.4 
to 0.42±0.06 at z ~ 0.2. It is more difficult to quantify the halo mass of objects in this sample 
because of the high false group detection rate (see §8.4). This sample is expected to contain 
a broad mix of subhalo masses, with very roughly half of the galaxies being in subhalos with 
masses of ~ 1 — 3 x IO^^Mq, and the other half, being contaminated by projections, residing 
in subhalos of smaller masses, similar to those in the non-group cluster galaxy sample. The 
more massive subhalos in this sample, along with the group cluster galaxy sample, could 
represent a total of 30 to 45% of all the the infalling cluster galaxies. They are likely the 
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main contributors to Butcher-Oemler effect seen in the outer radial bins. 

This rapid evolution of group populations at z ~ 0.3, combined with the red cores of 
clusters seen at higher redshifts and the bluer colors of galaxies in smaller groups, implies that 
galaxies in the most massive halos (i.e., in the cluster cores) have star formation quenched at 
an earlier time than those in moderately massive infalling groups. These, in turn, become red 
at an earlier time than the non-group galaxies residing in yet less massive subhalos (Figure 
12). This effect suggests that mechanisms responsible for the truncation of star formation in 
galaxies in more massive halos operate more efficiently, or become more effective in earlier 
epochs. 

Related to the more rapid evolution of fred in group galaxies is the 'pre-proc essing' effect 



on galaxies in infalling groups, suggested by number of previous studies (e.g., iBlake et al. 



2004J : iTanaka et al.ll2005l : ICooper et al.ll2008l : Mcintosh et al.l 120081 ). We clearly observe this 



in Figure [16], as the fred of galaxies in groups increases with decreasing redshift at all cluster- 
centric radii. The pre-processing effect is essentially a manifestation that the more massive 
subhalos in a cluster have a dominant effect on their galaxy populations, producing a strong 
Butcher-Oemler effect regardless of their distance from the cluster center. By z ~ 0.35, 
this 'processing' of galaxies in moderately massive subhalos results in galaxies which are 
redder than the non-group galaxies in similar local galaxy density regions located at similar 
cluster-centric radii. The fred of these galaxies, even at large tcl, becomes similar to that 
of the evolved galaxies in the cluster core by z ~ 0.2. While our sample demonstrates this 
effect in groups with subhalo masses larger than ~ 3 — 4 x 10^'^ Mq, it also likely occurs in 
subhalos 2 or 3 times less massive. This increase in the dominance of red galaxies in groups 
appears to be a significant, if not the primary, contributor of the Butcher-Oemler effect over 
the redshift range of 0.5 to 0.2. The increase of fred for the galaxies in both the non-group 
cluster galaxies and the cluster core galaxies is considerably smaller over this redshift range, 
in that the galaxies in the cluster core have already completed their evolution, and small 
subhalos (hosting the non-group cluster galaxies) outside the core still contain a substantial 
number of blue galaxies even ai z ^ 0.2. 



10. Conclusions 

We have investigated the environmental effects on the properties of the CNOCl cluster 
galaxies at 0.17 < z < 0.55 using both individual galaxies and those in galaxy groups, 
utilizing photometric redshift data covering a relatively large field out to 2 - 3 R2oo- We 
use local projected galaxy density, S5, and cluster-centric radius, rcL, in units of i?2oo as 
proxies for the local and global environment, and adopt the red galaxy fraction, fred, as the 
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measurement of galaxy population properties. We summarize our results as follows: 

1. ) The Butcher-Oemler effect is observed in all tcl and S5 bins between redshifts of 0.5 
and 0.2 (Figure [H]). However, the change of fred in the cluster core (less than 0.5 -R200) is 
minimal aX z < 0.4, suggesting that the truncation of star formation in the center of the 
most massive halos occurred at earlier epochs. 

2. ) We find that stellar-mass selected samples also show the Butcher-Oemler effect. Compar- 
ing the samples selected using different stellar-mass limits (at 1 -R200)) "we find the steepness 
of the Butcher-Oemler effect to be dependent on the stellar-mass limit. A sample selected 
using a higher stellar-mass limit has a milder change in fred compared to one with a lower 
stellar-mass limit. 

3. ) There are apparent correlations between S5 and fred in all our three redshift bins. 
However, the dependence of fred on local galaxy density becomes much weaker when the 
galaxy sample is separated by cluster-centric radius or group/non- group status, with the 
possible exception at the cluster outskirts. This suggests that local galaxy density may not 
be the primary environment parameter affecting the average galaxy population. 

4. ) We separate the cluster galaxies into three subsamples: cluster main-group, group cluster, 
and non-group cluster galaxies, roughly representing galaxies situated in the core of the 
massive dark matter halo of the cluster, in moderately massive subhalos, and in low-mass 
subhalos, respectively. These three samples have different fred values and rates of change 
of fred with rcdshift. Galaxies in the central massive cluster halo reach fred ~ 1 at a larger 
redshift than those in galaxy groups, which in turn reach fred ~ 1 at an earlier time than 
the non-group galaxies. This suggests that the dominant determinant of the epoch of the 
quenching of star formation is the mass of the dark matter subhalos in which the galaxy 
resides. The apparent correlation between E5 and fred is at least in part due to the correlation 
of high local galaxy density with the presence of a massive halo. The epoch when fred in a 
galaxy group reaches the evolved stage of fred ~ 1 depends on the halo mass. 

5. ) We find that galaxy groups falling into a cluster 'process' galaxies therein and turn 
them red at an earlier time; what is termed as 'pre-processing'. This is observed even after 
controlling for S5. Such 'pre-processing' is seen for groups at the outer tcl and the cluster 
outskirts. By z ~ 0.20, galaxies in groups have a high red fraction, similar to that in the 
cluster core, regardless of their cluster-centric radius. This more rapid evolution of galaxies 
in infalling moderately rich groups appears to be a significant contributor to the Butcher- 
Oemler effect. 

Some of our analyses and interpretations are tentative and suggestive, based on a small- 
size sample of 16 clusters and 89 cluster groups. The completion of the Red-sequence Cluster 



Survey 1 and 2 (RCS: lHsieh et al.ll2005l : lYee et al.ll2007l ) promises to provide similar data for 
thousands of clusters, providing greatly improved statistics to investigate the roles of dark 
matter halo environment on galaxy evolution. 
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Table 1: Properties of CNOCl Follow-up Data Set 
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Fig. 1. — The estimated photometric redshift for 2022 CNOCl galaxies with spectroscopic 
redshifts. The top two panels show the comparison between the spectroscopic and photomet- 
ric redshifts. The overall dispersion is a{zphot — ^spec) ~ 0.051. At z < 0.1, the photometric 
redshift measurements are affected by not having a sufficiently blue band. The bottom two 
panels plot Zphot — Zspec as functions of galaxy magnitude and color. The scatter is larger 
toward fainter magnitudes and blue colors: c^Zphot — Zgpec) ~ 0.048 and ~ 0.108 for Rc < 21.5 
and Rc > 21.5, respectively, and <j{zphot — Zgpec) ~ 0.041 for B — Rc> 1.8 and ~ 0.070 for 
B-Rc< 1.8 
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Fig. 2. — The photometric redshift distribution (open histogram) of all galaxies with 
rriR^ < 22.80 in the field of Abell 2390 field {z = 0.23). The background galaxy distri- 
bution estimated by the RCS CFHT patches (Hsieh et al. 2005) is shown as the dashed 
line. The cluster galaxies (hatched histogram), as determined by the group finding algo- 
rithm (see 2]), show a large excess in galaxy number in photometric redshift space at the 
cluster redshift. 
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Fig. 3. — The map of Abell 2390 with groups found using the pFoF algorithm. Galaxies 
selected in the cluster redshift space are plotted as dots. The symbols with different shapes 
denote members of cluster groups with A^^^^ > 8. The cluster main-group members are 
marked by open circles. The contours of the local projected galaxy density S5, calculated 
using on all cluster galaxies, are overlaid with levels of £5=20 and 40 gal/Mpc^. The dashed 
circle has a radius of 1 -R2oo- 
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Fig. 4. — The observed color-magnitude diagram of galaxies {M^^ < M^^ + 2.0) of galaxies 
i n the cluster main-group and galaxy groups in Abell 2390. The theoretical red-sequence 
( iKodama fc Arimotdll997l ). adjusted for calibration systematics (see §6.3), is overlaid in each 
panel as the solid line, while the dashed lines define the region used for red galaxies in the 
calculation of fred- Group galaxies in each panel are plotted using the symbol corresponding 
to the same group in Figure [31 Large symbols mark galaxies with M'^^ < M^^ + 1.5, which 



are used for the computation of fred- 
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Table 2: Exposure Times^ and Filters for Each Cluster 
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Fig. 5. — The red galaxy fractions for our 16 clusters using (1) galaxies in the identified 
cluster main-groups, displayed as solid circles; cluster galaxies (2) within 1 -R200! presented 
as open diamonds; and (3) at 1-1.5 -R20O) plotted as triangles. A linear fit between fred and 
cluster spectroscopic redshift is performed for each sample, with the slope and its uncertainty 
indicated on the plot. 
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Table 3: Properties of pFoF Determined CNOCl Clusters {DO^y = 0.25Mpc) 



Cluster R.A.(J2000) Dec.(J2000) ZpFoF,ci Nggi R^^^ fred 
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Fig. 6. — The estimated stellar mass M.,, function of M^'^ for Abell 2390 and MS1358+62. 
The dots are all cluster galaxies in the patch and the open circles mark those within I-R200 
of the cluster center. The vertical dashed lines indicate + 1.0 and M^^ + 1.5. The 
horizontal dashed lines mark log = 9.75, 10.2, and 10.6, respectively, increasing from the 
bottom. 
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Fig. 7. — The fred for galaxies within I-R200 of each cluster as a function of redshift. Results 
using three different stellar- mass (M^,) limits are shown. A linear fit between fred and clus- 
ter spectroscopic redshift is performed for each sample, with the slope and its uncertainty 
indicated on the plot. 
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Fig. 8. — The median stellar mass as a function of redshift for red and blue galaxies in the 
clusters. Galaxies are selected within I-R200 and in the absolute magnitude range + < 
M'^ < + 1.5. Red galaxy samples are represented by crosses, and blue galaxy samples 
by open squares. These median stellar-mass values provide relative estimates of the stellar- 
mass limits for the red and blue galaxies in the cluster sample at the luminosity limit of 
M^^ -|- 1.5. The horizontal dashed lines show the median of the values from all the clusters. 
The stellar-mass limits are similar for all the clusters in the sample. The small increases 
seen in the limits for the highest redshift clusters are due to slight incompleteness in absolute 
magnitude. 
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Fig. 9. — The fred as a function of tcl in three redshift bins, as indicated on the plot. The 
rcL values in each redshift bin are shifted slightly to provide clarity. Galaxies in all S5 
regions are included. 
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Fig. 10. — The fred as a function of S5 in three redshift bins, as indicated on the plot. The 
fred increases with increasing S5 in all three redshift bins. 
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Fig. 11. — The ' f red- z' trends in different environments, where vcl increases from the left to 
the right and S5 increases from the top to the bottom. The general trend of fred decreasing 
with increasing redshift is seen in all panels. 
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Fig. 12. — The 'fred-^b ' trends in different tcl and redshift bins. The trend of fred increasing 
with S5 is relatively weak compared to those in Figure [lOl This trend is most apparent at 
low redshift and at large tcl- 




Fig. 13. — The 'fred-fcL ' trends in different S5 and redshift bins. The trend with vcl is 
strongest in the ^5=0-10 bin and for the high-redshift subsample in the other S5 bins. 
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Fig. 14. — Local galaxy density S5 as a function of cluster-centric radius tcl in units of 
i?2oo using Abell 2390 cluster galaxies. The crosses, squares, triangles, stars, and diamonds 
represent group galaxies in different cluster groups (see Figure [3l). The thick solid circles 
are the 'non-group' cluster galaxies and those marked with large open circles are the cluster 
galaxies in the cluster main-group. The small dots are those in between our 'group' and 'non- 
group' categories. In general, the group galaxies populate regions of S5 relatively higher than 
those of the non-group galaxies. The solid curve is the projected NFW profile fitted using 
the cluster main-group and non-group galaxies. 
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Fig. 15. — The 'fred-z' trends for the cluster main-group, group cluster, and non-group 
cluster galaxies as indicated on the plot, using galaxies in all S5 bins. All group and non- 
group galaxies are located at tql < 3. The cluster main-group subsample has the largest 
fred, while the non-group galaxies exhibit a weak 'fred-z' trend. 
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Fig. 16. — The 'fred-f'CL ' trends in different redshift bins for the three subsamples indicated 
on the plot. Group cluster galaxies outside ~ I-R200 evolve much more rapidly than non- 
group cluster galaxies over this redshift range. 
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Fig. 17. — The 'fred-^5 ' trends in different redshift bins for the three subsamples indicated 
on the plot. 



